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The  high-quality  Fermi  LAT  observations  of  gamma-ray  pulsars  have  opened  a  new  window  to  understanding  the  generation 
mechanisms  of  high-energy  emission  from  these  systems.  The  high  statistics  allow  for  careful  modeling  of  the  light  curve  features 
as  well  as  for  phase  resolved  spectral  modeling.  We  modeled  the  LAT  light  curves  of  the  Vela  and  CTA  1  pulsars  with  simulated 
high-energy  light  curves  generated  from  geometrical  representations  of  the  outer  gap  and  slot  gap  emission  models,  within  the 
vacuum  retarded  dipole  and  force-free  fields.  A  Markov  Chain  Monte  Carlo  maximum  likelihood  method  was  used  to  explore 
the  phase  space  of  the  magnetic  inclination  angle,  viewing  angle,  maximum  emission  radius,  and  gap  width.  We  also  used  the 
measured  spectral  cutoff  energies  to  estimate  the  accelerating  parallel  electric  field  dependence  on  radius,  under  the  assumptions 
that  the  high-energy  emission  is  dominated  by  curvature  radiation  and  the  geometry  (radius  of  emission  and  minimum  radius  of 
curvature  of  the  magnetic  field  lines)  is  determined  by  the  best  fitting  light  curves  for  each  model.  We  find  that  light  curves  from 
the  vacuum  field  more  closely  match  the  observed  light  curves  and  multiwavelength  constraints,  and  that  the  calculated  parallel 
electric  field  can  place  additional  constraints  on  the  emission  geometry. 


1.  INTRODUCTION 

The  pulsar  emission  mechanism  is  not  well  understood. 
Magnetospheric  particle  acceleration  is  likely  responsible 
for  the  observed  emission,  but  the  emission  geometry  is 
unknown.  One  can  gain  some  insight  by  comparing  light 
curves  derived  from  geometrical  emission  models  with  ob¬ 
served  pulsar  light  curves.  Observations  of  pulsars  by 
the  Fermi  Gamma-ray  Space  Telescope  Large  Area  Tele¬ 
scope  (LAT)  [4]  have  shown  that  the  high-energy  emission 
likely  originates  in  the  outer  magnetosphere  [1],  We  sim¬ 
ulated  high-energy  light  curves  from  geometrical  versions 
of  two  standard  high-altitude  emission  models,  the  outer 
gap  (OG)  [13]  and  slot  gap  (SG)  [11],  and  from  a  third, 
modified  SG  model  with  azimuthal  asymmetry  in  emis- 
sivity  due  to  a  naturally  occurring  offset  dipole  (aSG)  [8], 
These  models  were  considered  within  two  field  geome¬ 
tries,  the  vacuum  retarded  dipole  (VRD)  and  force-free 
(FF)  [6]  fields.  We  compared  the  resulting  light  curves 
with  the  LAT  light  curves  of  the  Vela  pulsar  and  PSR 
J0007+7303,  the  CTA  1  pulsar,  to  constrain  the  systems’ 
geometries,  and  calculated  the  model-dependent  magni¬ 
tude  of  the  accelerating  electric  field  E\\  in  the  VRD  field. 


2.  LIGHT  CURVE  MODELING 

To  model  the  LAT  light  curves,  we  first  simulated  pul¬ 
sar  light  curves  from  geometrical  representations  of  the 
SG,  aSG,  and  OG  emission  zones  within  the  VRD  and  FF 
fields,  following  the  simulation  method  of  [7],  The  B  field 
defined  in  the  observer’s  frame  is  transformed  to  the  co¬ 
rotating  frame  (CF)  [5]  and  photons  are  emitted  tangent  to 
B  in  the  CF  prior  to  calculation  of  the  aberration.  We  as¬ 
sume  constant  emissivity  along  the  field  lines  in  the  CF. 
The  azimuthal  asymmetry  in  the  polar  cap  (PC)  angle  for 


the  aSG  model  is  calculated  for  each  inclination  angle  a 
as  in  [8]  (see  also  [9]).  For  a  given  a,  gap  width  w  (in 
units  of  open  volume  coordinates  rovc,  as  in  [7]),  and  max¬ 
imum  emission  altitude  r,  the  code  outputs  the  dimension¬ 
less  emission  intensity  and  the  minimum  and  maximum 
radii  of  curvature,  pmm  and  pmax,  emission  radii  rmjn  and 
rlnax,  and  local  field  magnitude  |/f|mln  and  |Z?|max,  at  all  ob¬ 
server  angles  7  and  rotation  phases  (ft. 

We  simulated  light  curves  for  a  fiducial  rotation  period 
of  0.1  s  on  a  4-dimensional  grid  of  a,  7,  w,  and  r.  Our 
simulation  resolutions  are  1°  in  a  for  the  VRD  field  and 
15°  for  the  FF  magnetosphere;  1°  in  7 ;  0.01  rovc  in  0  rovc  < 
w  <  0.3  rovc;  and  0.1  R\c  in  0.7  R\c  <  r  <  2.0  R\c,  where 
Ric  =  c/Q  is  the  light  cylinder  radius.  Emission  is  allowed 
out  to  a  cylindrical  radius  rcy i  =  0.98  R\c  for  the  OG  model 
and  rcyi  =  0.95  R\c  for  the  SG  models. 

LAT  light  curves  were  constructed  from  photons  in 
an  angular  radius  8  <  max[1.6  -  31og10(£’),  1.3]  from 
the  pulsar  [2].  The  PSR  J0007+7303  light  curve  has 
32  fixed-width  bins,  and  was  taken  from  [3],  The  Vela 
light  curve  has  140  fixed-count  bins  of  ~  3000  pho¬ 
tons  each.  The  background  of  PSR  J0007+7303,  195 
counts/bin,  was  found  using  the  Fermi  tool  gtsrcprob  as 
in  [3].  The  emission  in  the  off-peak  above  the  background 
level  was  assumed  to  be  magnetospheric  in  origin  (but 
see  [3]  for  details).  The  Vela  background  was  found  to 
be  204  counts/bin  using  the  off-peak  (phases  0.8-1,  where 
no  magnetospheric  emission  was  detected  in  our  spectral 
fits)  counts  in  the  energy-dependent  PSF  of  the  pulsar. 

We  used  a  Markov  Chain  Monte  Carlo  (MCMC)  maxi¬ 
mum  likelihood  routine  [14]  to  search  the  parameter  space 
for  the  combination  of  (a,  7 >  w,  r,  A <p)  that  best  repro¬ 
duced  the  LAT  light  curves.  The  fifth  parameter,  A (f>,  is  the 
amount  by  which  a  model  light  curve  must  shift  in  phase 
in  order  to  best  match  the  LAT  light  curve.  The  MCMC 
begins  at  a  random  point  in  parameter  space,  calculates  the 
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Figure  1:  Light  curve  modeling  results  for  Vela  and  PSR  J0007+7303.  The  absolute  best  fit  parameters  are  given  next  to  each  model; 
uncertainties  in  individual  parameters  vary,  but  are  typically  of  order  ~  10%.  In  all  panels,  the  light  curve  is  shown  in  black.  Blue 
curves  show  the  best  fit  OG  light  curve,  red  the  best  fit  SG  light  curve,  and  green  the  best  aSG  light  curve.  Vertical  dot-dashed  lines 
show  phase  zero  of  the  model  light  curve  of  corresponding  color;  the  phase  of  the  line  is  the  same  as  the  listed  A tp,  and  is  the  phase  of 
emission  coming  from  the  closest  magnetic  pole.  The  horizontal  purple  dashed  lines  show  the  background  count  level.  Panel  (a) 
shows  results  for  the  Vela  pulsar  with  the  VRD  held,  (b)  Vela  with  the  FF  held,  (c)  PSR  J0007+7303  with  the  VRD  held,  and  (d)  PSR 
J0007+7303  with  the  FF  held. 


likelihood  L,  and  moves  in  a  random  direction  to  a  new 
point  in  space  to  calculate  L  again.  If  L/Lformer  >  1  or 
>  rand[0, 1),  L  is  saved  in  a  chain;  the  routine  runs  until  the 
chain  contains  the  user-specified  number  of  steps.  Many 
chains  were  run  to  explore  the  whole  parameter  space.  We 
used  Wilks’  theorem,  AlnL  =  -Ay2/2,  to  calculate  L  at 
each  point  in  parameter  space.  To  perform  a  fit,  we  sub¬ 
tracted  the  background  level  from  the  LAT  light  curve,  re¬ 
binned  the  model  to  match  the  data,  and  normalized  the 
model  to  the  total  counts  in  the  LAT  light  curve.  For  each 
fit  in  the  vacuum  case,  we  ran  200  chains  with  20  steps 
each  to  adequately  sample  the  parameter  space.  For  the  FF 
fits,  we  ran  20  chains  of  20  steps  each  for  each  a.  Multi¬ 
wavelength  constraints  on  a  and  £  were  considered  after 
fitting.  Our  fit  results  are  shown  in  Figure  1.  Each  plot 
shows  three  instances  of  the  LAT  light  curve  with  the  best 
OG,  SG,  and  aSG  light  curve  superposed.  The  zero  phases 
of  each  model,  corresponding  to  the  nearest  magnetic  pole, 
are  given  by  the  vertical  lines,  and  the  background  levels 
with  horizontal  lines. 

We  find  that  for  the  (a)  VRD  and  (b)  FF  fields,  the  OG 
model  (blue)  statistically  fits  the  Vela  light  curve  better 
than  either  SG  model.  The  SG  (red)  produces  too  much 
off-peak  emission,  leading  to  high  x 2  values.  The  aSG 


(green)  reduces  the  background  significantly,  leading  to  a 
much  better  fit.  The  aSG  also  qualitatively  fits  the  peak 
emission  well,  as  its  main  peaks  are  the  correct  approx¬ 
imate  height  and  an  inner  peak  is  present.  All  three  fits 
have  0  close  to  the  value  determined  from  the  X-ray  torus 
geometry  [10],  £  ~  64°.  A(p  is  consistent  with  the  observed 
phase  lag  between  the  radio  and  y-ray  peaks  for  the  VRD 
models;  however,  for  the  FF  case,  A(/>  is  too  large. 

For  the  pulsar  in  CTA  1,  the  SG  models  fit  much  better 
than  the  OG  in  both  held  geometries.  There  is  no  con¬ 
straint  on  A (f>  due  to  the  lack  of  a  radio  detection.  The 
VRD  geometry  in  (c)  produces  much  better  fits  than  the 
FF  in  (d);  there  is  little  difference  between  the  SG  and  aSG 
due  to  the  small  a.  A  larger  FF  a  leads  to  a  larger  PC  offset, 
which  lowers  the  first  peak.  The  large  difference  between 
a  and  is  consistent  with  the  pulsar  being  radio-quiet  due 
to  geometry-the  radio  beam  would  not  cross  our  line  of 
sight  for  such  a  large  \a  -  £|. 


3.  CALCULATION  OF  En 

The  y-ray  spectra  of  pulsars  are  well  fit  by  an  exponen¬ 
tially  cut-off  power  law. 
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where  No  is  the  differential  flux.  Ft,  the  energy  scale,  I  the 
power  law  index,  and  Ec  the  cutoff  energy;  b  =  1  results 
in  a  simple  exponential  cutoff,  while  b  <  1  gives  a  sub¬ 
exponential  cutoff  and  b  >  la  super-exponential  cutoff. 
For  phase  averaged  spectra,  b  <  1  due  to  blending  of  Ec 
as  it  varies  with  phase  (e.g.  [2]),  while  b  is  consistent  with 
(and  is  fixed  to)  1  in  individual  phase  bins. 

In  current  models  of  pulsar  emission,  at  energies  above 
~  100  MeV  the  emission  is  dominated  by  curvature  radia¬ 
tion.  Particles  reach  the  radiation  reaction  limit  at  Lorentz 
factors  ycR  ~  107.  In  this  limit,  the  curvature  radiation 
cutoff  energy  Eqr  is  related  to  E\\  and  pc  by 


Ec  r  =  - — 7cr  =  0-32/Ic 
^  Pc 


(2) 


0.0  0.2  0.4  0.6  0.8 

Pulse  Phase 


Assuming  all  emission  with  E  >  100  MeV  is  due  to 
pure  curvature  radiation  and  that  the  VRD  is  the  true  B 
field  structure,  we  calculated  E\\  in  each  light  curve  phase 
bin.  We  used  the  simulated  minimum  radii  of  curvature 
(pc  =  pmin)  from  the  best  fit  VRD  light  curves  of  §3  and 
the  measured  cutoff  energies  (£cr  =  Ec)  in  each  phase  bin. 
The  cutoff  energies  are  given  in  [3]  for  PSR  10007+7303. 
We  updated  the  Vela  0.1-100  GeV  phase  resolved  spectral 
results  with  30  months  of  LAT  data,  following  the  method 
of  [2]  with  3000  pulsed  counts  per  bin,  and  used  our  mea¬ 
sured  cutoff  energies  for  the  Vela  E\\  calculation.  As  an 
example.  Figure  2  shows  the  measured  Ec,  the  simulated 
Pmin,  rm;n,  and  |B|max,  and  the  calculated  E\\  for  the  Vela 
pulsar  peak  emission  (phases  0  <  </>  <  0.8),  using  the  best 
fit  parameters  from  the  VRD  SG  geometry. 

We  explore  how  the  parallel  electric  field  varies  with 
emission  altitude.  We  have  calculated  E\\  in  each  phase 
bin  for  the  best  fit  vacuum  OG  and  SG  model  parameters. 
Because  the  value  of  |B|max  corresponds  to  rmin,  we  have 
plotted  £11  and  E\\l\B\max  with  minimum  emission  radius 
for  the  OG  and  SG  models  in  Figure  3. 

For  both  emission  models.  Vela  has  overall  a  ~  constant 
or  gradually  varying  E\\  with  altitude  (panel  (a)  of  Fig¬ 
ure  3),  which  is  expected  (e.g.  [2]).  In  panel  ( b ),  the  value 
of  E\\  is  compared  with  |B|max.  As  expected  for  an  E  field 
induced  by  B,  the  ratio  E\\/\B\max  <  1  for  all  rm jn  out  to 
the  light  cylinder  radius  (near  and  beyond  1  R\c,  the  vector 
components  of  B  are  less  certain  and  are  not  included  in 
this  calculation). 

For  PSR  J0007+7303,  the  magnitude  of  E\\  with  rlnm  is 
consistent  with  a  constant  (panel  (c)),  and  its  values  are 
similar  to  those  calculated  for  Vela.  Note  that  the  geomet¬ 
rical  parameters  obtained  from  the  light  curve  fits  are  very 
different  from  those  of  Vela,  leading  to  a  different  range  of 
r min  for  PSR  J0007+7303.  The  ratio  £||/|B|max  <  1  for  the 
best  fit  SG  model,  but  it  is  >  1  for  the  best  fit  parameters  of 


Figure  2:  Example  of  simulated  and  measured  values  described 
in  §3.  (a)  The  LAT  light  curve  of  Vela  (purple)  and  measured  Ec 
with  phase  (yellow  diamonds),  (b)  Simulated  pmin  (solid  blue) 
and  rmi„  (dashed  green)  for  the  best  VRD  SG  fit  parameters,  (c) 
Simulated  |B|max  (solid  cyan)  and  calculated  E\\  (red  squares)  for 
the  best  VRD  SG  fit  parameters. 

the  OG  model.  There  are  instances  in,  for  example,  a  non¬ 
ideal  magnetosphere  [10]  where  E\\  may  be  >  \B\.  In  the 
case  of  the  vacuum  field,  in  which  there  are  no  currents, 
the  only  source  of  £j|  is  induction  by  B.  and  therefore  E\\ 
cannot  be  larger  than  |_B|max.  In  this  particular  case,  then, 
we  find  that  the  combination  of  the  OG  model  and  VRD 
field  we  have  used  does  not  approximate  the  physical  envi¬ 
ronment  of  the  pulsar  magnetosphere  and/or  the  geometry 
of  the  emission  zone. 


4.  CONCLUSIONS 

We  have  evaluated  the  geometries  of  the  slot  gap  and 
outer  gap  emission  models,  and  the  vacuum  retarded 
dipole  and  force-free  magnetic  field  solutions,  by  compar¬ 
ing  the  simulated  light  curves  with  the  LAT  light  curves 
of  Vela  and  PSR  J0007+7303  and  finding  the  geometri¬ 
cal  model  parameters  that  best  reproduce  the  data  in  each 
case.  In  general,  the  OG  has  no  off-peak  emission  (this 
is  largely  responsible  for  the  OG  fitting  Vela  the  best  and 
PSR  J0007+7303  the  worst),  while  the  SG  models  do  a 
better  job  of  reproducing  wing  emission  but  over-predict 
the  off-peak  emission.  Introducing  azimuthal  asymmetry 
of  the  PC  angle  in  the  aSG  model  leads  to  a  reduction  in 
the  off-peak  emission  level,  improving  upon  the  SG  model 
light  curve  fits  within  the  VRD  field.  In  the  FF  field,  how¬ 
ever,  aSG  light  curves  tend  to  have  a  much  reduced  first 
peak,  leading  to  significantly  worse  light  curve  fits. 
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Figure  3:  Ey  (a)  and  the  ratio  of  £j|/|.B|max  (b)  for  the  OG  (blue  squares)  and  SG  (red  diamonds)  models  that  best  fit  the  Vela  LAT  light 
curve,  (c)  Same  as  (a)  for  PSR  J0007+7303.  (d)  Same  as  ( b )  for  PSR  J0007+7303. 


The  first  y-ray  peak  occurs  at  later  phase  in  FF  light 
curves,  so  values  of  A ft  are  larger  in  the  FF  magnetosphere 
than  in  the  VRD  field.  Physically,  Aft  cannot  be  larger 
than  the  phase  lag  between  the  radio  and  first  y-ray  peak 
unless  the  radio  beam  model  is  highly  contrived.  The  FF 
field  requires  a  Aft  larger  than  this  phase  lag  for  Vela.  This 
suggests  that  the  true  pulsar  magnetosphere  may  be  signif¬ 
icantly  different  from  the  FF  magnetosphere. 

For  Vela,  all  fits  within  the  VRD  field  have  ft  close  to 
the  expected  value  of  64°,  and  all  have  reasonable  values 
of  |a  -  4 1  such  that  the  radio  emission  is  observable.  Inter¬ 
estingly,  the  aSG  model  has  ft  closest  to  64°,  and  while  its 
X~  is  poor,  it  is  an  improvement  over  the  SG  and  qualita¬ 
tively  reproduces  well  the  major  features  (two  main  peaks 
and  inner  peak)  of  the  pulsed  emission.  The  FF  fits  also 
get  close  to  the  correct  if  value,  but  only  the  OG  has  an 
acceptable  |cy  -  ft |,  while  the  best  fit  SG  models  are  consis¬ 
tent  with  a  radio-quiet  pulsar.  Both  field  structures  lead  to 
large  | a  -  ft\  for  PSR  J0007+7303,  consistent  with  the  lack 
of  detected  radio  emission. 

We  calculated  the  model-dependent  £11  for  the  OG  and 
SG  geometries,  and  found  that  it  is  constant  or  slowly 
varying  with  emission  radius.  Comparing  £j|  with  |/J|rnax 
leads  to  the  interesting  result  that  for  PSR  J0007+7303, 
£jl  >  |B|max  for  the  OG  parameters  that  best  fit  the  LAT 
light  curve.  This  is  not  consistent  with  the  VRD  where 
\E\  <  |B|,  and  thus  disfavors  the  OG  model  in  the  VRD 
geometry  at  the  location  in  parameter  space  where  the  best 
light  curve  fit  is  found.  The  calculation  of  E y  can  therefore 
be  used  as  a  diagnostic  of  the  model  magnetosphere  and 
emission  geometry,  in  addition  to  the  light  curve  fitting. 

We  cannot  rule  in  favor  of  the  SG  or  OG  from  these  light 
curve  fits.  As  the  models  are  purely  geometrical,  we  expect 
to  reproduce  only  dominant  light  curve  features,  and  our 
statistical  fits  are  poor.  The  vacuum  field  produces  better 
light  curve  fits  than  the  force-free  field.  However,  we  note 
that  the  resolution  in  a  of  our  FF  models  is  much  lower 


than  in  the  VRD,  so  further  modeling  with  higher  resolu¬ 
tion  is  needed  to  confirm  this  result.  We  have  demonstrated 
that  light  curve  modeling  leads  to  constraints  on  the  geom¬ 
etry  of  individual  systems,  and  that  the  phase  lag  is  an  im¬ 
portant  diagnostic  in  comparing  magnetic  field  structures. 
The  model-dependent  calculation  of  E\\,  and  comparison 
with  the  model  |_B|max,  can  additionally  be  used  to  constrain 
both  the  magnetosphere  and  emission  geometry. 
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